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We evaluate the coupling constant for the lightest scalar glueball to decay to pseudoscalar meson pairs. The 
calculation is done in the valence approximation on a 16 3 x 24 lattice at (3 = 5.70 for two different values of 
pseudoscalar meson mass. 



A valence approximation calculation [jy using 
between 25000 and 30000 gauge configurations 
at each of several different (3 gives a value of 
1740 ± 71 MeV for the infinite volume contin- 
uum limit of the mass of the lightest scalar glue- 
ball. An independent calculation (|] with between 
1000 and 3000 configurations at several different 
/3, when interpreted || following Ref. Q, yields 
1625 ± 94 MeV for the lightest scalar glueball 
mass. Finally, an analysis in Ref. || suggests 
that the error in the valence approximation to 
the scalar glueball mass is probably less than 60 
MeV. The mass calculation with larger statistics 
then appears to favor /j(1710) as a scalar glue- 
ball, with /o(1590) a possible but less likely can- 
didate. A third alternative is that the scalar glue- 
ball is a linear combination of the two states with 
significant contributions from each. No other well 
established resonance consistent with scalar glue- 
ball quantum numbers lies within 200 MeV of the 
1740 MeV prediction. 

The interpretation of the scalar glueball as 
/j(1710) or some combination of /,/(1710) and 
/o(1590) is obscured, however, by two unresolved 
issues. First, if the scalar glueball decay width 
were large enough, this particle would not be eas- 
ily identified, in which case it need not be any well 
established resonance. Second, since the scalar 
glueball is a flavor singlet, one might expect its 
coupling constant to a pair of rj's to be the same 
as its coupling to a pair of K's or a pair of it's. In 
fact, the observed values of these three coupling 
constants are far from equal both in the case of 
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/j(1710) and in the case of / (1590). 

We are now in the process of calculating 
the scalar glueball's coupling to pairs of pseu- 
doscalars. Our preliminary results strongly sug- 
gest this particle has a width significantly below 
200 MeV. Thus we consider it quite likely that 
the scalar glueball has already been found in ex- 
periment. In addition, the coupling constants we 
obtain show a variation with the mass of the de- 
cay products consistent with the variation in both 
the /j(1710) and the / (1590) couplings. 

To evaluate the glueball decay couplings, con- 
sider a euclidean lattice gauge theory, on a lattice 
L 3 x T, with the plaquette action for the gauge 
field, the Wilson action for quarks, and the gauge 
field transformed to lattice Coulomb gauge. Let 
U[(x) for a space direction i — 1, 2, 3, and integer 
radius r > 1, be a smeared link field found by 
averaging the r 2 links running in direction i from 
the sites of the (r— l)x(r— 1) square oriented 
in the two positive space directions orthogonal to 
i starting at site x. Let V[Ax) be the trace of the 
product of U((y) and fj(y) around the boundary 
of an (r — 1) x (r — 1) square in the ij plane begin- 
ning at x. Let the zero-momentum scalar glueball 
operator g r (t) be the sum of the V[j(x) for all i,j 
and x in the time t lattice hyperplane. Define the 

quark and antiquark fields (x) and \E' r (x) to 
be Wilson quark and antiquark fields smeared by 
convoluting the local Wilson fields with a space 
direction gaussian with root-mean-square radius 
r 1 . Define the smeared pseudoscalar field 7r[ (x) 

with flavor index i to be \P (x)j 5 Xi^ r (x), where 
Xi is a Gell-Mann flavor matrix. Let frj (k,t) be 
the Fourier transform of 7r£ (x) on the time t lat- 
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tice hyperplane. We now fix r and r' to the values 
3 and \f§, respectively, used in the summary of 
results presented below, and omit the radius su- 
perscripts. 

At large time separations t, the vacuum ex- 
pectation values < 7f|(0, t)7Ti(0, 0) > and < 
7rl(k,t)TTi(k,0) > approach ri{L 3 exp(— E^it) and 
*g\L 3 exp{— E^ 2 t), respectively. Here k is any 
momentum with norm 2ttL~ 1 , E v \ is the mass 
of a pion, and E^ 2 is the energy of a pion 
with momentum 2irL~ 1 . In addition, in the va- 
lence approximation, < g(t)g(0)) > approaches 
rjgL 3 exp(—E g t), where E g is the mass of the light- 
est scalar glueball. At this point it is convenient 
to redefine g(t), TTi(0,t) and ni(k,t) by dividing 
each by the corresponding r\. 

Define the two pion zero momentum opera- 
tor fl(0, t) to be G- 1 / 2 ^ -^(0,0^(0,0). Let 
the two pion operator IT(fc, t) for nonzero k be 
6" 1 J2ik ^i(k,t)ni(—k, 0) where the sum is over 
the six ways of placing a vector k of length k 
in a positive or negative lattice coordinate direc- 
tion. Define the vacumm subtracted operators 

n s (fc,t) toben(M)- < n|ft(M)|n>. Let |i > 

be the lowest energy eigenstate in IT S (0, i) | O >, 
and let |2 > be the lowest energy eigenstate in 
n s (27rL _1 ,t)|fi >. The states |1 > and |2 > are 
both normalized to 1. Let E„„i and E nl72 be the 
energies of these two states, respectively. Define 
the amplitudes r)ij(t) 
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l- 3 <i\IL s {0,t)\n>, 
L~ 3 < i\U s (2nL- 1 ,t)\n >, 



(1) 

(2) 



For large t, T}%j{t) has the asymptotic form 
rjijexp^—E^jt). Connected three-point functions 
from which coupling constants can be extracted 
are now given by 



= <g(t g )U s (0,t 7r )>, (3) 
= <g{tg)f\ s {2nL-\t 7r )> . (4) 



If the quark mass, and thus the pion mass, is 
chosen so that E^-ni is equal to E g , the lightest 
intermediate state which can appear between the 
glueball and pions in a transfer matrix expression 
for Ti(t g> t„) is 1 1 > . Thus for large enough t g 
with t n fixed, Ti(< g ,<^) will be proportional to 



the coupling contant of a glueball to two pions at 
rest. If the quark mass is chosen so that E^i 
is equal to E g , however, the lightest intermediate 
state which can appear between the glueball and 
pions in a transfer matrix expression for T 2 (t g , t„) 
is still |1 >, not 1 2 >. This holds since 7712 (*) 
in Eq. (|J) has no reason to vanish. To obtain 
from T 2 (t g ,t 7r ) the coupling of a glueball to two 
pions with momentum 2irL~ 1 , the contribution 
to T2(t s ,i 7r ) arising from the |1 > intermediate 
state must be cancelled off. From the three-point 
functions defined in Eqs. (||) and (Q) we therefore 
define the amplitudes 

Sl(t g , < x ) = Ti(t g ,t w ) — — — TTT^itgi iff)) (5) 
?722 V"iT I 

Sa(* fl ,0 = T 2 {t g ,t„) - ^^(i,^). (6) 

»7ll(*7r) 

In S-zitgjt-n) the contribution of the undesirable 
|1 > intermediate state has been cancelled. In 
Si(t g ,tir) a contribution from the intermediate 
state 1 2 > has been cancelled. Although the sub- 
traction in Siitgjt^) is irrelevant for large enough 
t g , we expect that as a result of this subtraction 
Si (~tg> tir) w iU approach its large t g behavior more 
rapidly than does T\ (t g , t^) . 

An additional intermediate state which can also 
appear in a transfer matrix expression for either 
T^tgit-n) is the isosinglet scalar bound state of a 
quark and an antiquark. The relative contribu- 
tion of this state to either three-point function, 
however, is suppressed by the factor jj\t g — t^\, 
where [i is the coupling constant of this state to 
the glueball. Phenomenological arguments fol- 
lowing Ref. H strongly suggest that /i is less than 
60 MeV. In our main results below, removing the 
additional intermediate state would therefore give 
at most a 4% correction. 

At large t g and t^, the three-point functions 
become 



Si (tg, t^) — > 

V3Ai7/ii(l - r 12 r 2 i)L 3 



V^EgE^ 



1 V*g , ^7r) , 



S2\tg, t w ) 



3A 2 772 2 (1 - r 12 r 2 i)E : 



S 2 (tg,t 7T ). 



(7) 
(8) 
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Here Ai and A2 are the glueball coupling con- 
stants to pseudoscalar pairs at rest or with mo- 
mentum 2irL~ 1 , respectively. Up to a factor of 
— i these coupling constants are invariant decay 
amplitudes with the standard normalization con- 
vention. The factors rjij are given by the large t 
behavior of T)ij(t) as discussed earlier, and r%j is 
Vij/Vjj- The factors s^tg,^) are 

Si(tg,U) = ^ e Xp(-E g \t - tg\~ 
t 

E ni \t-U\-E„i\t\). (9) 

In Eq. (|^) we make the simplifying assumption 
that Emri is 2E 7t i. In our results to be presented 
below we actually use a slightly more complicated 
expression for Si{t g ,t v ) which takes into account 
small differences between E^i and 2E 7T i. 

To obtain values of A^ from Eqs. (f7j) and (jg) 
we need the amplitudes %(£). These we deter- 
mine from propagators for two pion states. Define 
the two pion operator II(ii , £2) from smeared pion 
fields, not Fourier transformed and not renormal- 
ized, to be 6 -1 / 2 ^ ""i(0, £1)^,(0, t 2 ). Define two 
pion propagators by 

C x {t 1: t 2 ) = 

<U{t 1 +2t 2 ,t 1 +t 2 )fl s (0,t 2 )>, (10) 
C 2 (h,t 2 ) = 

< n(ii + 2< 2 , t X + t 2 )fl s (2TT/L, t 2 ) > . (11) 

For large values of ti, these amplitudes approach 
Ci(ti,t 2 ) = 

Cnexp(-E n7Tl ti) + C 2l exp(-E^ 2 ti), (12) 
Cij = 

ViVijih^iiih) + v / 6r?2??y(i2)??i2(i2), (13) 

where the r\i are determined from the asymptotic 
behavior of single pion propagators as discussed 
earlier. 

At (3 of 5.7, we found that k of 0.1650 nearly 
gives EVtttI equal to E g , and k of 0.1675 nearly 
gives E^tt-2 equal to E g . For k of 0.1650 and t from 
to 5, we determined rjij(t) with Eqs. Jl2] ) and 
( p^| ) applied to Ci{t\,t 2 ) from an ensemble of 100 
independent gauge configurations on a 16 3 x 40 
lattice. For k of 0.1675 and t of to 5, we deter- 
mined T)ij(t) from Ci(ti,t 2 ) using 870 indepen- 
dent gauge configuratons. The T)ij(t) turn out to 
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Figure 1. Decay couplings. 

be nearly independent of t. In each case 7700 is 
statistically consistent with 1.0, yyn is between 
1.0 and 1.1 and rji 2 and — 7721 are positive and 
less than 0.1. 

We then evaluated Ai for k of 0.1650 and A2 
for k of 0.1675 using Eqs. (0) and (||) applied 
to Si(tg,t„) and S 2 (t g ,t 7r ) found with an ensem- 
ble of 7200 independent gauge configurations on 
a lattice 16 3 x 24. We obtained statistically sig- 
nificant results for tg — t^ of 0, 1 and 2 and t n of 1, 
2, 3, 4, and 5. The Ai are nearly constant, within 
statistical errors, across these 15 points. Figure [l] 
shows predicted coupling constants fitted to the 
data with t g — t % of 1 in comparision to observed 
decay couplings for decays of /j(1710) to pairs of 
rj's, K's and n's. The horizontal axis gives the 
squared mass of the decay products. Masses and 
decay constants are shown in units of the p mass. 
If couplings to r/'s, K's and n's are determined 
from our predicted couplings by linear interpola- 
tion in the squared mass of the decay products, we 
predict a total width for glueball decay to pseu- 
doscalar pairs of 52^| MeV, in comparison to 
84 ±23 MeV for /./(1710). 

We would like to thank Steve Sharpe for finding 
an error in an earlier version of this work. 
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